Tunneling into Multiwalled Carbon Nanotubes: Coulomb Blockade and Fano 

Resonance 
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Tunneling spectroscopy measurements of single tunnel junctions formed between multiwalled car- 
bon nanotubes (MWNTs) and a normal metal are reported. Intrinsic Coulomb interactions in the 
MWNTs give rise to a strong zero-bias suppression of a tunneling density of states (TDOS) that 
can be fitted numerically to the environmental quantum-fluctuation (EQF) theory. An asymmetric 
conductance anomaly near zero bias is found at low temperatures and interpreted as Fano resonance 
in the strong tunneling regime. 
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Coulomb blockade (CB) has been studied intensely in 
a multi-junction configuration, in which electron tunnel 
rates from the environment to a capacitively isolated "is- 
land" are blocked by the e-e interaction if the thermal 
fluctuation is below the charging energy Ec = e^/2C and 
the quantum fluctuation is suppressed with sufficiently 
large tunnel resistance Rt » Rq = /2e. In the case 
of a single-junction circuit, the understanding of CB is 
less straightforward. The Coulomb gap, supposedly less 
signiflcant for the case of low-impedance environments, 
should be established only if the environmental impen- 
dence exceeds Rq. 

In single- walled carbon nanotubes (SWNTs), the re- 
duced geometry gives rise to strong e-e interaction. In- 
deed, CB oscillations and evidence of Luttinger liquid 
(LL) have been observed 0. In contrast to SWNTs, in 
which only two conductance channels are available for 
current transport, MWNTs with diameter in the range 
of (i=20-40 nm have several tens of conductance chan- 
nels. The energy separation of the quantized subbands, 
given by Ai? = hvf/d, is about 13-26 meV taking the 
Fermi velocity Vf — 8 x 10^ ni/s. This value is about 
an order of magnitude smaller than that of SWNTs. Ex- 
periments indicate that MWNTs are considerably hole- 
doped, thus a large number of subbands, on the order 
of ten, are occupied. Observations of weak localization 
(3, electron phase interference effect and universal 
conductance fluctuations 3] support the view that low 
frequency conductance in MWNTs is contributed mostly 
by the outmost graphene shell and is characterized by 2D 
diffusive transport. In addition to the phase interference 
effects, a strong e-e interaction has also been observed 
in MWNTs. Pronounced zero-bias suppression of the 
TDOS has been observed several times in the tunneling 
measurements 4, 5, 6]. Moreover, the TDOS shows a 
power-law, i.e., v{E) ~ E", which resembles the case 
of a LL. It is noteworthy that in the EQF theory, for a 
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single tunnel junction coupled to high-impedance trans- 
mission lines, such a scaling behavior is also predicted at 
the limit of many parallel transmission modes 0. The 
physical origin of these power laws is the linear disper- 
sion of bosonic excitations that are characteristic both 
for LL, which is a strictly ID ballistic conductor, and a 
single tunnel junction connected to a 3D disordered con- 
ductor. In the latter case, the quasiparticle tunneling is 
suppressed at 1^ <C e/2C, therefore the charge is trans- 
ported with ID plasmon modes. The fact that MWNTs 
have many conductance modes, together with the obser- 
vation of a crossover from power-law to Ohmic behavior 
at higher voltages @, suggests that the EQF theory is 
more appropriate to describe the observed TDOS renor- 
malization. However, most of these measurements were 
done in multi-junction conflgurations. A single tunnel 
junction measurement is needed to further clarify this 
issue. 

In this Letter, millimeter-long CVD-grown MWNTs 9 
are measured by a cross-junction method. The MWNTs 
samples are composed of loosely entangled nanotubes of 
diameters between 20 to 40 nm that are roughly parallel 
to each other and up to 2 mm long. The single tunnel 
junctions are formed by crossing a very thin (< 1/im) 
MWNT bundle with a narrow strip of metal wire fabri- 
cated on an insulating substrate. We explored different 
electrode materials including Au, Cu, Sn, and Al. In the 
case of Sn and Al, a small magnetic field is applied to 
suppress the superconducting state below Tc. No obvi- 
ous change of device characteristics attributable to the 
choice of metals is observed. With such a cross-junction 
configuration the measured conductance is contributed 
exclusively by the tunnel junction, since the current is 
passed along one arm of the MWNT/metal and the volt- 
age is measured along the other, non-current carrying 
arm. Thus the device can be understood as a small num- 
ber of single junctions in parallel. Despite the simplicity 
of their fabrication, we find that the devices are very 
stable and sustain several cooling cycles without appar- 
ent change of characteristics. More than 20 samples are 
measured, all yielding strong zero-bias suppression of the 
TDOS. 
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FIG. 1: (a) dl /dV as a function of V in dimensionless units 
measured at T = 0.12, 0.25, 0.78, 1.2, 1.7, 4.5, 10, 20 K 
(sample 990530s6, curves are offset for clarity). Dashed lines 
are the fits of EQF theory, (b) DC current simultaneously 
measured at 0.12 K. Dashed line is the numerical fit. Dotted 
line shows a Coulomb offset of e/2C (c) The temperature 
dependence of the exponent a for three samples (o: 990530s6; 
A: 990316; V: 990320). 



The current I{V) and G = dl/dV are calculated by 
a golden-rule approach incorporating the environmen- 
tal influence by P{E), the probability for a tunneling 
electron to lose energy to the environment. P{E) can 
be calculated from its Fourier transform, the phase cor- 
relation function J{t). In the transmission-line model, 
J{t) is determined by the total environmental impedance 
Zt{ijj) = [iLuC + Z~^{ijj))~^ , where Z{lli) is the external 
environmental impedance. We apply the method in Ref. 
to evaluate P{E) from an integral equation without 
the need of going to the time domain. The I-V character- 
istics are then calculated for finite temperatures and arbi- 
trary junction impedance. The parameters that appear 
in the numerical calculation are the damping strength 
a = Z(0)/Rq, the inverse relaxation time lurci f^nd the 
quality factor Q = lorc I^s with the resonance frequency 
of the undamped circuit los = {LC)^^^'^ . It is noted that 
the quality factor Q plays only a minor role as an addi- 
tional adjustable parameter. Therefore it is not impor- 
tant to the fit, and is always set to unit. 

We then try to fit our experimental data with the EQF 
theory. Since Ec and Rt are determined by the high- 
voltage data, the isolation resistance Riso = Z{Q) is the 
only adjustable parameter in our fitting. In contrast to 
the case of Ref. jlfl] , where Riso is formed by ideal Ohmic 
and temperature-independent resistors, Riso in our case 
are provided by the resistive impedance of the MWNTs 
themselves, which should be temperature-dependent. In- 
deed, from the fitting, we see an evolution of a from 0.1 
at 20 K to 0.25 at 1-4 K and then a saturates (Fig. Ic). 



TABLE I: A partial list of characteristics of the samples. 



Sample 


Rt (k!^) 


Ec (eV) 


C (aF) 


Rrso (kO) 


a 


990316 


18.4 


0.019 


4.2 


3.3" 


0.25" 


990320 


48.3 


0.014 


5.7 


3.3" 


0.25" 


990530s6 


8.6 


0.01 


8.0 


3.3" 


0.25" 


990530sl 


10.0 


0.02 


4.0 


4.6" 


0.35" 


990202 


9.1 


0.004 


20.0 


3.5" 


0.27" 



"value acquired at T = 1.2-4.3 K. 



Note that dimensionless units are used with G normal- 
ized by l/i?t and voltage normalized by e /2C. Therefore 
the number of the single tunnel junctions has no effect 
on the fitting, and our measurements of different sam- 
ples can be directly compared. We find that for differ- 
ent samples with scattered characteristics (see Table I), 
at low temperatures the exponent a reaches a univer- 
sal value of 0.25-0.35, which agrees with previous multi- 
junction measurements EH- It yields R iso — 3.3-4.6 
kfi, which is roughly a constant for different samples. 
This coincidence is not accidental but reflects the in- 
trinsic electrodynamic modes of the MWNTs, which can 
be modelled as an ideal resistive LC transmission line: 
R = (i'/C)^/^ with the kinetic inductance estimated 
as L' = RQ/2Nvf « 1 nH/^m for iV « 10 - 20 modes 
and the capacitance C' w 20-30 aF//j,m. Therefore, the 
"environment" with respect to the single-junctions in our 
devices is provided by MWNTs themselves, not by the 
external circuits. 

As mentioned previously, for an isolated single tun- 
nel junction with many parallel transport modes, 
the EQF theory predicts a power-law asymptotics: 
G(F,T)/G(0,T) = fiV/T) - |r[i(a + 2) + 
2^]/r[i(a + 2)]r[l -f where G(0,T) is the 

zero-bias conductance, and T{x) is the Gamma func- 
tion 0. For eV/2-KkBT 1, a voltage power law 
G(y, T) ~ V" is expected. Note that in the above scal- 
ing function, exponent a is the only adjustable param- 
eter. In Fig. 2, we can observe such scaling behavior. 
Here f{V/T), represented by soHd line, is calculated tak- 
ing a = 0.27. The inset shows the original data in which 
the dashed line is calculated numerically using the same 
exponent a = 0.27. 

CB in MWNTs is recently described with a micro- 
scopic theory in Ref. In the 2D diffusive regime, 
the exponent for tunneling into the bulk of a MWNT is: 
a = {R/ hv[)D)ln{l+V{)U[)) . Here R is the tube radius, Uq 
is the intra-tubc Coulomb interaction, D — Ujt/2 is the 
charge diffusivity and the "bare" DOS vq = N/Ahvf with 
N Ri 20. If we use u/t « 60 nm from the magnetoresis- 
tance data 12], then under the condition oi Uo/hvf ^ 1 
the above formula yields a « 0.25, which agrees with our 
experiment. 

Besides the voltage power law phenomenon, which is 
characteristic of strong Coulomb interactions, another 
energy scale - the discrete energy levels due to electron 
confinement - emerges as temperature decreases, and it 
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FIG. 2: Scaled conductance G{V, T) /G(0, T) of another sam- 
ple (990202). The inset shows the original dl/dVvs.V with 
the dashed line calculated with exponent a = 0.27. 



modifies the tunneling spectra. Unlike the case of a two- 
junction configuration, where the electrons are confined 
by the two contacts if the nanotube is clean enough, in 
a single-junction configuration the electrons can be con- 
fined by disorders when the MWNTs are "dirty" , such 
that the impedance of the local environment of the junc- 
tion is larger than Rq . A stacking mismatch between ad- 
jacent walls and other structural imperfections are possi- 
ble sources of disorders in MWNTs, resulting in discrete 
energy levels. 

By cooling down the devices to below 1 K, we observe 
that G develops a narrow resonance-like anomaly at very 
low bias (Fig. 3). The asymmetric anomaly builds up 
consistently as temperature decreases and even shows a 
dip structure. The line shape resembles that of a Fano 
resonance. 

Unlike the Coulomb structure, which is caused by 
static e-e interaction, a Fano resonance arises from an 
e-e exchange interaction between two interfering scatter- 
ing channels: a discrete energy level and a continuum 
band. It has recently been rediscovered in mesoscopic 
systems such as semiconductor quantum dots, SWNTs, 
etc. 0, 0, 0, . Taking a quantum dot as an ex- 
ample, it can be considered as a gate-confined droplet 
of electrons with localized states. The coupling of the 
dot to the leads can be tuned to control the system to 
enter different transport regimes: If the dot is weakly 
coupled with the environment, a well-established CB de- 
velops. The charge transport is suppressed except for 
narrow resonances at charge degeneracy points. When 
the tunnel barriers become more transparent, the dot 
enters the Kondo regime - i.e., below a characteristic 
Kondo temperature Tk, spin-fiip co-tunneling events in- 
troduce a narrow symmetric TDOS peak at Ep that can 
be interpreted as a discrete level. If the coupling is strong 
enough, the interference between this discrete level and 
the conduction continuum gives rise to an asymmetric 
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FIG. 3: Asymmetric resonance features of dl /dV seen in two 
samples: 990530sl (a) and 990530s6 (b). Dashed lines are the 
fits of Fano's formula. 



TABLE II: Characteristics of the two samples in Figure 3. 



Sample 


T (K) 


<? 


£o (meV) 


7 (meV) 


Tk (K) 


990530s6 


0.25 


0.53 


0.43 


0.32 


1.86 


990530s6 


0.12 


1.75 


0.57 


0.29 


1.68 


990530sl 


0.25 


0.92 


0.3 


0.65 


3.77 


990530sl 


0.08 


1.97 


0.43 


0.56 


3.25 



Fano resonance. In Ref. such a crossover from 

a well-established CB through Kondo regime to Fano 
regime has been clearly observed. 

We find that the asymmetric resonance curve of G can 
be fitted by the Fano's formula: G ~ (e qf/{e^ + 1). 
Here q is the so-called asymmetry parameter, e = {eV — 
eo)/(7/2) is the dimensionless detuning from resonance, 
and 7 is the FWHM of the resonance. As expected, the 
asymmetric parameter 5 as a measure of the degree of 
coupling between the discrete state and the continuum 
increases when the temperature drops. 

The Kondo temperature, estimated from the relation 
7 = 2kBTK, is consistent for each device at different tem- 
peratures, and agrees with the observed FWHM of the 
resonance. If Kondo physics truly exists, then G at volt- 
ages near the resonance should show non-monotonic tem- 
perature dependence around Tx 15]. Indeed, the G — V 
curves in Fig. 3 exhibit such behavior: G at the peak po- 
sition first drops with T and then rises up below ~ 1 K 
(plotted in Fig. 4b). Moreover, we find that a perpendic- 
ular magnetic field gives rise to effects that are two-fold 
(Fig. 4a and 4c): First, the background conductance in- 
creases monotonically with the applied B. Second, the 
dip seen at zero field gradually disappears and the reso- 
nance turns into a nearly symmetric peak at high field, 
similar to the pattern observed in semiconductor quan- 
tum dots 111. 

The magnetic field effect can be explained as follows: 
First, adding a flux should change the amplitude and/or 
the phase for the resonant channels and therefore break 
down the coherent backscattering and increase the for- 
ward transmission through the channel. Second, the 
magnetic field can destroy the interference between the 



4 




T (K) V (mV) 



FIG. 4: (a) The dependence of zero-bias dl /dV on perpen- 
dicular magnetic field, (b) The temperature dependence of 
dl/dV measured at different voltages, (c) Effect of perpen- 
dicular magnetic field. From bottom to top: B = 0,1,3,6 
Tesla (sample 990530s6). 

resonant and nonresonant paths, transforming a resonant 
dip into a peak. 

Since Kondo physics is historically interpreted as the 
interplay between the d-orbitals of magnetic impurities 
and the conduction continuum, we have to ehminate the 
possibihty that the observed Fano resonance comes from 



residual traces of the Fe/Si catalyst in the MWNT sam- 
ples. It has been shown that magnetic impurities in 
MWNTs cause an enhancement of thermoelectric power 
h^, which is absent in our careful TEP measurements 
[iSj. Furthermore, no Fe signature can be detected in the 
body of the MWNT bundles within the instrumental res- 
olution in the energy dispersion X-ray and TEM studies. 
Therefore the observed Fano resonance must be an inher- 
ent property of MWNTs, which can be virtually treated 
as quantum dots strongly coupled to the leads. We no- 
ticed that recently Buitelaar et al. have also seen clear 
traces of CB and Kondo resonance in MWNT SET de- 
vices Their observed Kondo temperature Tk = 1.2 
K is in good agreement with our results. 

In summary, strong zero-bias suppression of the TDOS 
is observed in MWNTs' single tunnel junctions that can 
be explained well by the EQF theory. The observed ex- 
ponent a « 0.25-0.35 is found to be consistent for all the 
samples. Similar to the case of an open quantum dot, in 
low-impedance junctions we find that a Fano-resonance- 
like asymmetric conductance anomaly builts up below 
mV energy scales. It seems that MWNTs provide us a 
good laboratory to study the interplay between strong 
e-e interaction and disorder scattering. 
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